Two different quasiparticle scattering rates in vortex line liquid phase of layered 

d-wave superconductors 
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We carry out a quantum mechanical analysis of the behavior of nodal quasiparticles in the vortex 
line liquid phase of planar d-wave superconductors. Applying a novel path integral technique we 
calculate a number of experimentally relevant observables and demonstrate that in the low-field 
regime the quasiparticle scattering rates deduced from photoemission and thermal transport data 
can be markedly different from that extracted from tunneling, specific heat, superfluid stiffness or 
spin-lattice relaxation time. 



In recent years, the physics of nodal quasiparticles 
in planar d-wave superconductors such as the high-T c 
cuprates has attracted a lot of both, theoretical and ex- 
perimental, attention. The spectroscopic and transport 
properties of these Dirac-like quasiparticles with a linear 
dispersion have been studied in quite some detail, and 
elaborate analyses of various mechanisms of elastic scat- 
tering in the uniform superconducting state have been 
carried out, including the effects of potential [1], Kondo- 
like [2] , and extended impurities [3] as well as twin bound- 
aries [4]. 

In the mixed state, the recent quantum mechanical 
generalization of the earlier semiclassical approach [5] 
proposed in Rcf. [6] allows one to account for the non- 
uniformity of the local d-wave order parameter A p (r) = 
A cos(26* p ) exp(i(/)(r)) by means of a singular gauge trans- 
formation from the physical electrons c CT (r) of spin a to 
the new fermionic quasiparticles. Unlike electrons, the 
latter are subject to the effective magnetic field with zero 
mean which, besides the physical field, also includes the 
supercurrent circulating outside vortex cores. 

Applying the gauge transformation of Ref. [6] to the 
electronic states with energies small compared to the 
maximum gap one can represent them in terms of the 
Nambu operators creating the auxiliary fermions with 
the momenta near the nodes of A p (r) 



c CT (r) 
d CT (r) 



E ■ 

n=l,2;± 



(1) 



where n = 1,2 labels the pairs of the opposite nodes, 
while the choice of the phases </>A,s(r) is only restricted 
by the condition 4>a(t) + <Ab(i") = <K r )- 

In the case of a regular vortex lattice, the representa- 
tion (1) was used to demonstrate that the structure of 
the quasiparticle energy spectrum is that of the energy 
bands, rather than the Landau levels [7]. 

In the present paper, we extend the analysis based 
on the representation (1) to the experimentally well- 
documented vortex line liquid (VLL) phase [8] where 
the vortices are distributed totally randomly due to their 
strong pinning by columnar or other defects. 

Any disorder, including that induced by random vor- 
tices, is expected to affect the behavior of the d-wave 



quasiparticles most strongly at the lowest energies, pos- 
sibly resulting in the complete quasiparticle localization 
which, however, is still awaiting for experimental confir- 
mation [9]. In what follows, we focus on the ballistic 
regime of quasiparticle energies large compared to the 
localization scale (see below) which is readily accessible 
by a number of standard probes such as angular-resolved 
photoemission (ARPES), thermal transport, tunneling, 
specific heat, muon spin rotation (/iSR), and spin- lattice 
relaxation. 

In the quantum-mechanical approach of Ref. [6], the 
combined effect of external magnetic field H = V x A(r) 
and swirling supercurrent characterized by the super- 
fluid velocity v s (r) = |(V<Aa + Ws) — f A does not 
amount solely to the semiclassical Doppler shift e — > 
e — ki?v s (r) of the quasiparticle energies [5]. The latter 
is to be complemented by the vector potential, a(r) = 
Tj(V(j)A — V0b), that couples to the quasiparticles via 
their momentum, k — > k — a(r), and accounts for the 
quantum mechanical Berry phase corresponding to their 
Bohm-Aharonov (BA) scattering by the vortices. Thus, 
the complete Hamiltonian of the non-interacting nodal 
quasiparticles contains both the scalar- and the vector- 
like random terms 



,(r) + v^iijpi - cii(r))]Vv- (2) 



In Eq.(2), we used the 4x4 representation for the ma- 
trices 7 M = ((72, i&i, 10^3) <8> 03 acting in the space of 
the Dirac bi-spinors composed of the Nambu spinors: 

In order to keep our discussion and formulas relatively 
simple we consider the case of isotropic quasiparticle dis- 
persion and use the units where Vi = hc/2e = ks = 1. 
Moreover, because of the predominantly small-angle na- 
ture of the quasiparticle scattering by the vortices we 
choose to neglect the processes of inter-node scattering. 
While anticipating that neither of these simplifying as- 
sumptions will affect our main conclusions, we recognize 
that a combination of the above factors in the case of the 
real cuprates may give rise to an additional one-to-two 
dimensional crossover regime [10]. 
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As in the previous studies of the VLL phase [11], we 
average over different vortex configurations by assuming 
the Gaussian distributions 



(g(»,o)(<0g(„,a)(- < 0) = ™Ka)(q) (Sij - , 



(3) 



for both g„ = v s and g a = a. The striking difference 
between w v (q) = a/(q 2 + a) and w a (q) = a/q 2 which 
are both proportional to the areal density of vortices 
a = 2ttH reflects the presence of screening for the scalar- 
like Doppler potential and its absence for the vector-like 
BA scattering. The longitudinal (oc <}Wj/q 2 ) part of the 
correlator (3) proves to contribute negligibly to all the 
quantities of interest (besides, it is suppressed by the 
Coulomb interactions). 

In contrast to the previous analyses which focused 
solely on the effect of Doppler scattering [11], we find 
that vortex disorder has a profound effect on the quasi- 
particle spectrum which can not be adequately modelled 
by a constant quasiparticle width. In order to illustrate 
this point, we apply the standard self-consistent Born 
equation 



S(e,p) 



f e 

J (27r) 2 q 2 - 



e + S(e,q) 



( £ + S(6,q)) 2 



W(<,,a)(P + q) (4) 



to the separate contributions of the two scattering mecha- 
nisms towards the total quasiparticle width ImE t ,+Im£ a . 

In the case of scalar disorder Eq.(4) yields Im£„ oc a 1 / 2 
for small energies and momenta (e,p < a 1 / 2 ) while at 
max(e,p) 3> a 1 / 2 it behaves as oc a/max(e,p) which 
dominates over the scattering by the vortex cores whose 
rate is estimated as oc a/ A. 

At first sight, the effect of the BA scattering may seem 
to be much stronger, since a naive solution S a of Eq.(4) 
with the singular kernel w a (q) is plagued with a logarith- 
mic infrared divergence of the momentum integral. 

In order to avoid this spurious divergence which stems 
from the non-gauge invariant nature of the auxiliary 
fcrmion propagator one has to proceed directly with com- 
puting the manifestly gauge-invariant (retarded) Green 
function of the physical electrons 



= E - 

2;± 



p oo 

G R (e,r)= dte iet (c(t,r)J(0,O)} = Y 

J ° „=!,: 

(V>„(r) expH / (v s + aff 3 ® l)dr']^ n (0)). 
Jc 



(5) 



It turns out that the exponential decay of this function 
(see below) makes it largely independent of the contour 
C which can then be chosen as the straight path between 
the end points r and 0. 

To compute the amplitude (5) we apply the path- 
integral representation of Ref. [12] to the propagator of 
the auxiliary Dirac fermions. First, for a fixed vortex 



configuration, we cast Eq.(5) in the form of a functional 
integral over the space-time coordinate r(r) and the con- 
jugate momentum p(r) parameterized by the proper time 



oo f-r{r) — r 

s '--i" ^ - / a~ I DrDpe iSoM 



S K (e,r|v s ,a)= / dr 

JO Jr(0)=0 



exp[i^ dr'(k F v s + ^a) -i J (v s + a)dr')], (6) 

where we dropped, for the sake of compactness, the sum 
over the nodal points and introduced the free fermion 
action 



S [r] = jf dT '[e7o + p(J^-7)]- 



(7) 



Averaging over the disorder variables v s and a with 
the use of Eq.(3) results in the electron Green func- 
tion G R (e, r) = (Q R (e, r|v s , a)) which is given by Eq.(6) 
where, instead of the exponential phase factor, the inte- 
grand contains a product of two attenuation factors 

W M [r(r)} =exp[-- j ^ J q ^ J q ^ 

u J (r 1 )u l (r 2 ) W ; ( ^ Q) (q) e ^ r ^)- r ^»], (8) 

with u v = vf — r/r and u a = dr/dr — r/r. Thus, the 
presence of the exponent of the line interal taken along 
the contour C in Eq. (5) strongly reduces the effect of both 
a and v s , as compared to the case of the gauge- variant 
propagator of the auxiliary fermions. This observation 
seems to have been overlooked in the earlier studies of 
this and related problems where the phase factor in ques- 
tion would either not appear at all [13] or be averaged 
separately from the fermion propagator computed in a 
different approximation [14]. 

Proceeding along the lines of the previous analyses of 
the problem of non-relativistic fermions subject to a ran- 
dom vector potential [15], one can show that in the bal- 
listic regime (which in the present case is defined by the 
condition e 3> a 1 / 2 ) the path integral (6) is dominated by 
the fermion trajectories which only slightly depart from 
the straight line r (r) = vr. Evaluating the factor W v 
for such a trajectory, one obtains 



W w [r(r)] w exp 



, (9) 



while the integral in W a turns out to be proportional to 
the so-called Amperian area of the closed contour com- 
posed of a fermion trajectory r(r) and the "return" path 
— ro(r). Although this purely geometrical term vanishes 
for the saddle-point trajectory r(r) = i"o( T )> its expan- 
sion to first order in the transverse deviation r± (r) yields 
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W [r(r)]«cxp[-| jf dr'\r ± (r')\}. (10) 

By analogy with the non-relativistic problem studied in 
Refs. [15] the path integral (6) with the Wi factors given 
by Eqs.(9) and (10) can be related to the resolvent of the 
Schroedinger equation describing the transverse motion 
of the Dirac fcrmion 



d 2 + (e 2 - q 2 ) + {ae/qf {\x\ + x Q ) f + iasignx 
xg(e, q\x, x') = 5(x — x'), 



(11) 



where x = (y/W$ - VP) 2 /2a 1 / 2 . 

By analogy with the results of Refs. [15] the averaged 
physical electron propagator G R (e, p) can be obtained by 
convoluting the kernel l/(p 2 — q 2 ) 3 ^ 2 with the solution 
of (11) taken at x = x' = which is given by the formula 



<?(£,#, 0) = 



-i -l 



Tx^ u ^ 



(12) 



x=0 



where U± — U (a±, X[xo ± x]) is the parabolic cylinder 
function of the parameter a± — (e 2 — q 2 ± ia)/X 2 and 
A= (2ia\e/q\) 1 / 2 . 

Turning now to the applications of Eq.(5), we first dis- 
cuss the electron spectral function satisfying the disper- 
sion relation G R (e, p) = / A(e', p)de' /n(e- e' + iS). Near 
the maximum, |e 2 — p 2 | < a, it can be shown to take the 
form 



A(e,p) 



(e% - P7)v^ 



[ (e 2_ p 2 ) 2 + / g2 (1 + ( | c/ p|_ 1) 4 /4)] 3/4 



(13) 



where (3 ~ ira/8, thereby demonstrating a replacement 
of the bare pole by a branch cut of the function (z — z ) 3 ^ 2 
resulting from the above convolution procedure. Accord- 
ing to Eq.(13), the decay of the electron propagator in 
real space (G R (e,r) oc e(r jr) 1 l 2 e~ Vr for r ^> 1/r) is gov- 
erned by T(e) oc a/e which should be thought of as the 
actual (energy-dependent) quasiparticle width. 

This direct experimenta prediction can be tested by 
performing ARPES measurements in the VLL phase of 
the cupratcs under the weak-field conditions (VH <C 
T <c A). 

Also, comparing Eq.(13) with the estimate for Im£„ 
obtained from Eq.(4), we conclude that in the ballistic 
regime both the BA and the Doppler scattering mech- 
anisms appear to be equally important, contrary to the 
conclusions drawn in Ref. [16]. 

Next, we compute thermal conductivity given by the 
averaged product of two electron propagators 



(e/T) 2 de 



cosh 2 (e/2T)7 2tt 



dr 



Tr< 7l ^(e,r) 7l e*(e,-r)>. 

(14) 



The corresponding path integral reads as 

poo /■ri I 2(n,2)=±r 

(jg A (e,r)W R (e,-r)) = / d Tl dr 2 

JO Jr 1>2 (0)=0 

II Dr a Dp^e i§ °^e iS °M ^ Wl (r a -rp), (15) 

a.f3— 1,2 i—a,v 

where the factors Wi(r a — rp) with a ^ (3 account 
for the vertex corrections, alongside the self-energy ones 
(a = 0). Thus, the path- integral method of comput- 
ing the Dirac fermion conductivity is capable of pro- 
ceeding beyond the conventional (non-crossing and fan- 
shaped) ladder series of the vertex corrections to the bare 
fcrmion bubble (the latter suffice only if the number N 
of the Dirac species is large, while in the d-wave problem 

^ = £ CT i = 2). 

Upon integrating over the "center of mass" variables 
ri + V2 and pi + P2 and rescaling the ones describing 
relative motion, we again arrive at Eq.(ll). This time 
around, it is formulated in terms of the transverse rela- 
tive coordinate (r\± — r2_i_), and its solution yields the 
average {G A g R ) cx e(r 4r /r) 1 / 2 e - r "' r , where T tr = 2T is 
now playing the role of the momentum relaxation rate, in 
agreement with the above estimate for the quasiparticle 
spectral width. Plugging this asymptote into Eq.(14) we 
find that in the low-field regime the thermal conductivity 
behaves as 



15 H' 



(16) 



While being in agreement with the estimate based on the 
kinetic equation k oc Tv/T tr proportional to the linear 
density of states (DOS) v cx T and T tr oc H/T, Eq.(16) 
is strikingly different from the result (k oc T 2 / if 1 / 2 ) that 
one would obtain by naively assuming that the quasi- 
particle width remains constant (r' oc H 1 / 2 ) up to the 
energies e^T (cf. with Refs. [11]). 

The analysis of the data of Ref. [17] taken in 
Y Ba-zCuzOs.gg shows that Eq.(16) should be expected 
to hold for T < 30K and 0.1 < VW/T < 1 where the 
vortex-induced T tr dominates over the other mechanisms 
of scattering, including potential impurities. 

Our approach also enables one to compute other av- 
erages, such as (G R g R ) oc e(T tr /r)^ 2 e 2ler - Tt - r , which 
controls the effect of vortex disorder on superfluid stiff- 
ness measured by fiSR 



p s (0)-p s (T) 



irk, 



-In, / detanh(^) (17) 



/ 



drTr(£*( e , !•)£%, -r)) 



21n2, 



-T + 



H , T 2 



Notably, in contrast to the spectral and transport charac- 
teristics whose behavior is determined by the structure of 
the electron spectral function near its maximum, Eq.(17) 
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is governed by the overall momentum integral of the so- 
lution of the two-particle analog of Eq.(ll). 

This behavior is common amongst the thermodynamic 
quantities associated with the averages of local bi-linear 
operators VcrMVv ( r ) which are invariant under the 
gauge transformation (1). In fact, such averages as, e.g., 
G R (e, 0) must be computed differently, since now the tra- 
jectories contributing to the path integral (6) may devi- 
ate very strongly from the semiclassical one, r (r) = 0, 
although the mean square of the distance over which a 
typical trajectory r(r) ventures from the origin still scales 
quadratically with time, (r 2 (r)) oc t 2 . For 1/A < r this 
allows one to evaluate the damping factors as 



Wj[r(r)] w exp 



T 



Plugging (18) instead of (8) into (6) and computing the 
resulting (quadratic) path-integral, we arrive at the cor- 
rection to the linear DOS corresponding to the clean limit 



V {e)=TmTr[%G*{e,0)]=^?(± 



-H, (19) 



where a 2 H oc Hln(A 2 /H) and F{x) = tt 1 / 2 x Erf (a;) + 
cxp(— x 2 ). The effect of vortex disorder is most pro- 
nounced at small energies e < H 1 / 2 , and it appears to 
be stronger than in the semiclassical (Doppler-only) ap- 
proximation (cf. with [18]). 

Directly, this DOS correction can be extracted from the 
tunneling conductance G(V) oc v{V). Indirectly, it can 
also be manifested through the correction to electronic 
specific heat 



C(T) 



(e/T) 2 v(e)de 18C(3)T 2 



' H 



o cosli (e/2T) 



IGttT 2 



(20) 



Notably, the correction AC oc H 2 ln( A 2 /H)/T 2 is 
smaller than the result (AC oc H) obtained under in 
the situation where the inter-vortex repulsion is stronger 
than random pinning, and therefore the VLL is partially 
ordered [18]. 

Although the smallness of the disorder-induced term 
in Eq.(20) might hinder its detection, an alternate possi- 
bility is offered by the spin-lattice relaxation time 



Ti(T) K , 



v 2 (e)de 



cosh 2 (e/2T) 3 



'H 



3tt 5 / 2 : 



(21) 



where we dropped the overall prefactor proportional to 
the ion-specific matrix elements. 

In summary, we carried out a fully quantum mcchnical 
analysis of the quasiparticle properties of the VLL phase 
of layered ci-wave superconductors. We demonstrated 
that both the semiclassical Doppler shift and the intrin- 
sically quantum mechanical BA scattering have compa- 
rable effects on all the observables. Our path-integral 



approach enabled us to identify the energy-dependent 
effective quasiparticle width T(e) oc H/e describing the 
near-maximum (e 2 w p 2 ) behavior of the distinctly non- 
Lorentzian electron spectral function which can be di- 
rectly measured by ARPES. 

Also, we exposed the striking difference between the 
two rates: F tr ~ r(T) oc H/T displayed by the transport 
characteristics and T' ~ T(i/ 1 / 2 ) oc H 1 / 2 manifested by 
tunneling and thermodynamic quantities. 

In contrast, the real-space averaging procedure applied 
in the previous studies of the disordered vortex states 
[18] is bound to deliver the latter rate, because, focusing 
solely on the energy distribution, it does not faithfully 
represent the momentum dispersion of the averaged elec- 
tron spectral function. 

We emphasize that the origin (that is, a non-trivial en- 
ergy dependence) of the different apparent quasiparticle 
rates, as revealed by the different measurements, must 
be distinguished from the conventional juxtaposition of 
quasiparticle lifetime versus transport time in, e.g., pho- 
toemission and transport experiments. 

We expect that, albeit derived under a number of sim- 
plifying assumptions, our main conclusions are robust 
against including such factors as spatial anisotropy and 
non-linearity of the fcrmion dispersion as well as intcr- 
node scattering and, most importantly, modifying the 
distributions (3) with the purpose of describing partially 
ordered "vortex glass" states [18]. 

Lastly, the results of this paper can also be used to 
describe the effects of thermal phase fluctuations con- 
trolled by a oc T in the pseudogap phase of the cuprates 
[13,14]. To this end, wc predict that different mea- 
surements may return different values of the effective 
quasiparticle width, which might explain the inconsis- 
tency between the widths deduced from the tunneling 
and ARPES data [13]. 

This research was supported by the NSF under Grant 
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